Electromagnetic properties of graphene junctions 
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A resonant chiral tunneling (CT) across a graphene junction (GJ) induced by an external elec- 
tromagnetic field (EF) is studied. Modulation of the electron and hole wavefunction phases ip by 
the external EF during the CT processes strongly impacts the CT directional diagram. Therefore 
the a.c. transport characteristics of GJs depend on the EF polarization and frequency considerably. 
The GJ shows great promises for various nanoelectronic applications working in the THz diapason. 
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Unconventional electromagnetic properties of 
monoatomic layer graphene [J Q show remarkable 
perspectives for scientific research and novel electronic 
applications. An external electromagnetic field (EF) af- 
fects the electric charge transport [H across the graphene 
junction (GJ) substantially. The mechanism of such a 
strong influence implies an intrinsic origin, since the 
external EF modulates the wavefunction phase of the 
charge carriers during their transmission across the GJ 
shown in Fig. 1(a). The charge carriers in graphene 
behave as relativistic massless 'chiral fermions' (CF) 
characterized by a 'pseudospin'P, The pseudospin 
relates P, Q the CF particles to a certain sublattice A 
or B and is responsible for various unconventional d.c. 
electronic and magnetic properties of the monoatomic 
layer graphene Q. Significance of the pseudospin is 
proclaimed in the graphene junction (GJ) formed by a 
graphene sheet with two electrodes attached as shown 
in the main Fig. 1. The electrodes can either be 
metallic or made of the graphene itself. The potential 
barrier is induced by the gate voltage Vq from a Si 
gate. The chiral tunneling (CT) strongly depends on 
the azimuthal angle (j) between the two-dimensional 
electron momentum p = [hk, hq) and x axis (see Fig. 1). 
The CT is intrinsically related to a relativistic nature 
of CF since the electron and hole wavefunctions are 
interconnected in the monoatomic graphene [E 0. At 
some selected angles = 0„ the CT probability becomes 
ideal (T = 1). Such a resonant tunneling occurs due to a 
constructive quantum interference between an incoming 
electron (e) and a hole {h) moving inside the graphene 
barrier in a reverse direction and is characterized by 
the same pseudospin as sketched in Fig. 1 (Klein 
paradox^, Q). The interference pattern resulted from 
the chiral tunneling (CT) is very sensitive to the phase 
difference ip between e and h wavefunctions. Since an 
external electromagnetic field (EF) affects the phase ip 
directly, a thorough study of that quantum interference 
mechanism between EF, e and h is necessary. 

This Letter addresses the EF-induced directional reso- 
nant chiral tunneling across the GJ. We expect that angu- 
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FIG. 1: The graphene junction (GJ) controlled by the gate 
voltage Vg being exposed to the external electromagnetic field 
(EF). Two different barrier configurations (a) with additional 
dielectric sublayers and (b) without them. Below we show the 
corresponding graphene barrier profiles. 



lar dependence of the differential tunneling conductivity 
a {(j)) {4> is the angle between the electric current j and the 
x-direction indicated in Fig. 1) is strongly influenced by 
the EF since the latter modulates the wavefunction e — h 
phaseshift ip inside GJ. In this way the a.c. fleld adjusts 
the resonant tunneling condition, causing a directional 
shift of the resonances. Initially we compute the steady 
state CT (0,Vsd) curves for a GJ with either metallic or 
graphene electrodes biased by the source drain voltage 
VsD- The results obtained for the steady state are then 
utilized for studying of the a.c. properties on the next 
stage. In particular we will see that the EF suppresses 
the CT current in the straightforward direction {(p — 0) 
due to an angular redistribution of the d.c. current. The 
d.c. current increases to its resonant value at finite angles 
4> ^ instead. Below we focus on a weak a.c. field limit 
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and discuss origin of the intrinsic noise. Here we examine 
two different setups of the voltage biased graphene bar- 
rier. One setup (quoted as the rectangular barrier) cor- 
responds to a graphene barrier stripe 2 separated from 
the electrodes 1 and 3 [see Fig. 1(a)] by two thin di- 
electric sublayers at a; = and x ~ D. Since resistivity 
of the dielectric sublayer is several orders of magnitude 
higher that of the graphene stripe, the bias voltage Vsd 
drops mostly on the formers [see the profile below Fig. 
1(c)]. For the above reasons the d.c. electric field £ is 
negligibly small inside the graphene section of the bar- 
rier. In this case the CF wavefunction ^ inside the bar- 
rier graphene stripe is well approximated by plane waves. 
Another setup constitutes a bare monoatomic graphene 
stripe which contacts the attached graphene electrodes 
immediately as sketched in Fig. 1(d). The source drain 
bias voltage Vsd m this configuration drops entirely on 
the middle graphene stripe [see its profile below Fig. 
1(d)], which corresponds to a trapezoidal shape of the 
graphene barrier tilted proportionally to Vsd- For the 
GJ having finite dimensions, the motion of CF is quan- 
tized. The quantization imposes additional constrains on 
the directional tunneling diagram. 

The Dirac equation for fermions is written as 

- ihv {{dx ®i)dx + [oy ® T^) dy) * + eC/ (a;) = 

(1) 

where v = c/300 is the massless fermion speed, e is the 
CF excitation energy, di and fk are the Pauli matrices, 
{i. A:} — 1 ... 3, the barrier potential U {x) is induced by 
the gate voltage Vg. In the steady state, when the EF is 
off, the electric current is fully suppressed when Vsd < t/o 
(for typical gate voltage Vg = 1 V and the Si02 thick- 
ness d = 300 nm one finds[6] Uq ^ 2 meV). The CF 
excitation energy in the leads (x < and x > L) reads 

£p = eVoo ± hvy^k^ + = hv {kprn ± ^/k'^ + ■ Set- 
ting £p = for a normal incidence (g„ = 0) one writes 
eVoo ± hvkpm = or eVoo ~ fivkprn- If the leads are 
metallic, one writes k — \J {{sp — fi) /hv)^ — q^. The en- 
ergy inside the graphene stripe is £p — eVc^hv \J k^ + q^^. 
Inside the graphene barrier the CF wave vector is fc = 

\J {{sp — eC/o) /fiv)'^ — qn- The transverse quantization 
conforms to the boundary conditions at y = and 
y = W. The transmission amplitude tf, is conveniently 
computed using a simple approach, which represents the 
CF wavefunctions ^ as plane waves for a rectangular bar- 
rier and via Airy functions for a trapezoidal barrier. In 
particular, the scattering state for a trapezoidal barrier 
is constructed as 

* = 6'(-a;)[xn,feAi(fc, x) + r„x«,-fcBi(A:, x)] 
+e{x - i)t„Xn,fcAi(fc, X- L) + e{x)e{L - x) 

X [a„x„,fcAi(fc, a;)e'^ + /3„x„,_fcBi(fc, x)e^^] (2) 

where Ai(fc,a;) and Bi(fc,a;) are the Airy functions, 
and we also introduced the auxiliary functions 




FIG. 2: The steady state transparency |t(£,i?l')|^ versus < 
e < 2.5 and the incidence angle < <j) < n/2 for the rectan- 
gular and trapezoidal barriers. The contour plot of ctsd [in 
units of e^vpN (0)] versus Vsd (in units of the barrier height 
eUo) (a) for a rectangular barrier, and (b) for a trapezoidal 
barrier. 



Xn,fc(y) = a„|T) ® (IT) + ^n,fc|i))e"?"^ +aUi) ® 
{zn,k\]) + \l))e'^-y +5„|T) ® {zn,k\]) + \i))e-'^-y 
+K U> ® (IT) + ZnM li)) e"*«"2' where (g) is the Kronecker 
product, IT)^ = (1 0) and = ( 1 ) are 1x2 

matrices, Zn,k = ±{k + iqn) I \/k'^ + 9^, where ± signs 
apply to conductive (valence) bands, k is positive for 
conductance band and negative for the valence band, the 
factor Zn^k satisfies the identity Zn,kZn,-k = —1- Permit- 
ted values of the angle x inside the graphene stripe of a 
finite width W are obtained from boundary conditions 
along the y-direction, py = q„ = nn/W, W being the 
graphene stripe width. This gives the quantized angle Xn 
inside the graphene stripe as Xn = arctan [nn/ (kpW)]. 

The electric current between the electrodes 1 and 3 [see 
Fig. 1 (c,d)] is computed implementing methods Pzl, Isl, 1^. 
as j = 2TrveN (0) / de [Gf (e, p) - Gf (e, p)] where 
p = {hk,hq), V — 8.1 • 10^ m/s is the Fermi velocity 
in graphene, N (0) = 3.8 nm^^eV"^ is the electron 
density of states at the Fermi level, {e, p) = 

Ep IV I' e^«"''e''=^ {2np I) S {e - Ep + Sr,3eVsiy) 
is the Keldysh Green function r is the elec- 
trode index, Up is the distribution function of 
electrons with a 2D momentum p = {hk,hqn) 
assuming q is quantized. Then one obtains 
J - in/2)evN (0) / de Ep l^e |' ((2np - 1) J (e - + eV) 

- {2np -l)5{e - Ep)) = nevN (0) E„ / dE |ij,„|^ («£,„ - ns-ev,n) 
where N (0) is the electron density of states at the 
Fermi level. The steady state transparency diagrams 
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T = |t (e, 0)1 versus < e < 2.5 and the incidence angle 
< (/) < 7r/2 for a rectangular and a trapezoidal barrier 
are shown in Fig. 2. In the contour plots Fig. 2(a,b) one 
may distinguish two types of the quantization, (a) A 
conventional quantization in the lateral direction and (b) 
a quantization taking place in the straightforward direc- 
tion. The latter quantization is highly unconventional 
since it exists for the chiral ferniions only. Corresponding 
scattering states are sharply angular-dependent, which 
is evident from the same contour plot. The 3D plots 
in Fig. 2(c,d) indicate numerous resonances of ideal 
transparency [dark areas and spots visible in Fig. 2(a,b) 
corresponding to T — 1] occurring below the barrier 
when VsD < 2 (in units of the barrier height Uo) at finite 
incidence angles 0^0. 

If an external a.c. field is applied to the graphene junc- 
tion, its influence strongly depends on the field intensity 
and polarization. The fermion wavefunction phase differ- 
ence across the junction becomes modulated by the a.c. 
field ip (x, t) ^ Tp {x, t) exp {ieVi^c /* dt cos fit j . Below we 
discuss various a.c. field polarizations for the rectangular 
and trapezoidal barrier configurations. First we address 
the d.c. biased rectangular barrier shown in Fig. 1(c). 
If the electric field vector E of an external a.c. field is 
polarized as E = {E^, 0, 0), one may account for the a.c. 
field influence by setting Vsd Vsd + eVac cos fir, where 
Kc = E^D, E^ = -Vx'^ - A^,, $ and being the 
scalar and vector potentials. That approximation gives 
^ (x, t) ~* ^ (x, t) ■ expjieKc /* dr cos fir}. In this con- 
figuration the a.c. field does not affect the barrier shape, 
it simply splits the quantized scattering resonant levels 
only, i.e. Eq Eq ± nfi where Eq is a resonance en- 
ergy, n is the number of photons. That is the case when 
EF modulates the quantum mechanical phase difference 

across the junction. Using the equality llj e*"*''"^* = 
X]„ (a) e™^', one obtains the photon- assisted tun- 
neling (PAT) differential conductivity as cr (Vsd) ~ 
Y.n Jri {eVac/hfl) ctq (Vsd " nHfl/e) where ctq (Vsd) is the 
steady state differential tunneling conductivity computed 
above. Vac is amplitude of the a.c. bias voltage in- 
duced across the GJ by EF. The photon-assisted tunnel- 
ing is spectacularly proclaimed for the quantized scatter- 
ing states. A qualitatively different phenomena emerges 
if the external a.c. field applied to GJ modulates the 
barrier's height and shape. The barrier transparency T 
in that case is immediately affected by the a.c. field, po- 
larized, e.g., as E = {0,0, Ez). For a weak field Vac « 
Uq, the PAT correction to the d.c. tunneling current 
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{2n, 



where the angular-dependent transmission amplitude 
*e^n obtained from non-stationary boundary conditions 
at x = and x = L for the two configurations of the d.c. 
biased graphene barrier shown in Figs l(c,d). Since the 
EF modulates phase of the transmitted CFs, it modifies 
also the condition of the resonance chiral tunneling. The 
a.c. field modulates the phase difference between the 
incident electron (e) and the reverse-propagating hole 
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FIG. 3: The EF-induced excess angular-dependent trans- 
parency T\ of the rectangular (a) and trapezoidal (b) 
graphene barriers versus Vsd. (c) Fano factor -F(Vsd) of the 
intrinsic Poisson noise for three different lengths L — 5, 15, 30 
of the trapezoidal barrier, (d) The EF-induced deflection of 
the electrons and holes during the chiral tunneling. Here Ti 
and F are expressed in dimensionless units while Vsd is in 
units of the barrier height eUo- 



(h) waves propagating in opposite directions deflecting 
them inside the barrier [as sketched in Fig. 3(d)]. This 
causes splitting of the CF momentum inside the barrier 

as hk± — ±\J {{e ± ftfi — eU{))/hvY — sin^(0) where 
typical frequencies are in the THz diapason. The wave- 
function phaseshift inside the barrier is determined as 



(p± = arctan 



esin 4>/{hk±) 



where 6 is the electron in- 



cidence angle in the electrode attached to the graphene 
stripe. That means that the directional diagram strongly 
depends on the external field frequency fi (typically 0.5- 
50 THz) and magnitude Vac- In this way one gets a 
unique electromagnetic sensing devise, which tunning re- 
quires no energy excitation of the charge carriers, unlikely 
to most of the known sensors nowadays. 

The angular-dependent excess transmission Ti(e, <j)) in- 
duced by an a.c. field, which modulates the graphene 
barrier height Uq is shown in Fig. 3 where the plots (a) 
and (b) correspond to the rectangular and trapezoidal 
barriers. The a.c. field tends to suppress tunneling in 
the longitudinal direction [cj) = 0) while creates new res- 
onances with T = 1 at finite angles (jmeqO. 

Intrinsic noise in the monolayer graphene junction 
originates as follows. The thermal noise produced by 
phonons emitted in the electron-phonon collisions is de- 
pleted. The depletion is governed by quantum mechan- 
ical selection rules for the chirality conservation. That 
maintains the d.c. transport being ballistic up to room 
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temperatures. The intrinsic mechanism of the thermal 
noise depletion is as follows . Matrix element of the 
electron-phonon collisions according to Ref. [H, is 
Mppi oc {p \M {x)\p') cos (0pp' /2) where (j)pp' is the angle 
between the initial and final states. The phase factor 
cos {(f>pp' /2) ensures suppression of the electron-phonon 
and electron-impurity collisions. Therefore the charge 
transport remains ballistic up to room temperatures. An- 
other intrinsic noise (Poisson noise) arises due to the 'Zit- 
terbewegung' effect caused by a jittering motion of the 
change carriers when electrons are randomly converted 
to holes forth and back. The intrinsic noise character- 
istic are illustrated by the Fano factor i^(VsD). The 
Fano factor is computed using the above expressions for 
the junction transmission coefficient tpp' and the electric 
current. The noise persists even in the zero tempera- 
ture limit. Such a Poisson noise is characterized by the 
Fano factor F = J2nTn{l- T„) / J2n Tn, where r„ is the 
transmission probability in the n-th channel and the sum- 
mation is performed over all the conducting channels (in 
our setup this means just integration over 0). From the 
plot F(Vsd) shown in Fig. 3(c) for different length of the 
graphene barrier L = 5,15,30 [in units of vh/ {2TreUo), 
which is typically 1-6 nm] one infers that the Poisson 
noise has a sharp resonance for long junctions {L > 30) 



and substantially decreases when Vsd ^ Uq. 

In conclusion we computed the electric current across a 
monoatomic layer graphene junction in conditions when 
an external electromagnetic field is applied. We find that 
the a.c. properties of GJ strongly depend on the EF 
frequency and amplitude. The EF induces an angular 
redistribution of the tunneling current. Such a deflec- 
tion comes rather from a quantum mechanical phase shift 
than from an inelastic excitation by the EF. That en- 
sures no heating is involved during the absorption. The 
a.c. current induced by EF across the GJ has a sharp an- 
gular dependence, which potentially can be exploited in 
sensor nanodevices of the external electromagnetic field. 
An experimental observation of such a spectacular EF- 
induced deflection would provide strong supporting evi- 
dences for the chiral fermion concept in the monoatomic 
graphene. The EF-induced CT resonances improve the 
device switching speed and enhance the frequency res- 
olution. The phenomena considered above have a great 
potential for building of fast-switching transistors and EF 
sensors working in the THz diapason. 
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